Why air bubbles in water glow so easily 



Michael P. Brenner 1 , Sascha Hilgcnfcldt 2 , and DetlefLohse 2 

1 Department of Mathematics, 

Massachusetts Institute of Technology, 
Cambridge, MA 02139, USA 



2 



Fachbereich Physik, Philipps-Universitat Marburg, 
Renthof 6, D-35032 Marburg, Germany 



Abstract. Sound driven gas bubbles in water can emit light pulses (sonolumi- 
nescence). Experiments show a strong dependence on the type of gas dissolved in 
water. Air is found to be one of the most friendly gases towards this phenomenon. 
Recently, Lohse et al. (1996) have suggested a chemical mechanism to account for 
the strong dependence on the gas mixture: the dissociation of nitrogen at high tem- 
peratures and its subsequent chemical reactions to highly water soluble gases such 
as NO, NO2, and/or NH3. Here, we analyze the consequences of the theory and 
offer detailed comparison with the experimental data of Putterman's UCLA group. 
We can quantitatively account for heretofore unexplained results. In particular, we 
understand why the argon percentage in air is so essential for the observation of 
stable SL. 



1 Introduction 

Gaitan et al. (1990) discovered that a gas bubble levitated in water by a strong 
periodically modulated acoustic field can emit bursts of visible light so intense 
as to be observable to the naked eye, a phenomenon called single bubble 
sonoluminescence (SL). Previous workers had observed sonoluminescence in 
multibubble cavitation clouds: when a liquid is subjected to large negative 
pressures, the liquid rips apart into a cloud of bubbles, which then collapse 
violently. Frenzel and Schultes (1934) placed a photographic plate in the 
vicinity of a cavitating bubble cloud, and discovered that it was exposed to 
a low level of radiation. Single bubble SL is distinguished from multibubble 
sonoluminescene not only by the much higher light intensity, but also by the 
fact that the bubble remains in the oscillating state for billions of oscillation 
periods (days) without dissolving or changing its average size. 

Soon after Gaitan's initial discovery, Putterman's research group at UCLA 
began to uncover a myriad of interesting dynamical properties of single bub- 
ble SL (Barber and Putterman (1991), Hiller, Putterman, and Barber (1992), 
Hiller et al. (1994), Hiller and Putterman (1995), Barber et al. (1994), Bar- 
ber et al. (1995), Lofstedt, Barber, and Putterman (1993), Lofstedt et al. 
(1995), Weninger et al. (1995)). The width of the light pulse is less than 
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50 picoseconds. The accuracy of the flashes differs by less than 50 picosec- 
onds from cycle to cycle. The phenomenon shows strong dependence on the 
experimental parameters as the forcing pressure amplitude P a , the mass con- 
centration °f dissolved gas in the liquid far away from the bubble, and the 
temperature of the liquid. SL is only observed for P a ~ 1.1 atm - 1.5 atm (for 
an ambient pressure of P = 1 atm) for all gases. The dependence on is 
more subtle. For pure argon gas stable SL is observed only in a small window 
of extremely low concentration around c^ r /c ~ 0.4%, whereas for air the 
water has to be degassed only down to c^/co ~ 10 — 20% to obtain stable 
SL. Here, Co is the saturation concentration of the gas in the liquid. Unsta- 
ble SL, characterized by a "dancing bubble" and by oscillations in the phase 
and intensity of the light pulse, is observed in both cases (air and argon) for 
larger Coo/co- For pure nitrogen bubbles no stable SL is observed at all, only 
very weak unstable SL. Recently, Weninger et al. (1995) also observed single 
bubble SL in nonaqueous fluids (alcohols) where the phenomenon shows a 
strong dependence on the type of liquid. 

The major challenge in understanding the SL experiments is to determine 
how much of the phenomenology is hydrodynamic, how much is connected 
with the atomic properties of the gases, and how much stems from chemistry. 
To allow for detailed comparison between experiment and theory Hilgcnfeldt, 
Lohse, and Brenner (1996) calculated phase diagrams which depict domains 
of different bubble behaviors. They must depend on parameters which are 
easily experimentally adjustable. These are the forcing pressure P a and the 
gas concentration c x . As emphasized first by Hiller et al. (1994) the ambient 
radius i?o (i-e., the radius of the bubble under standard ambient conditions; 
i?o ~ 5/im in the SL experiments) is not an experimentally controllable 
parameter but adjusts itself dynamically. 

We will see that much of above phenomenology can already be understood 
from the hydrodynamics and the stability of the bubble. Three types of in- 
stabilities will turn out to be important: (i) Shape instabilities of the bubble 
(Plesset (1954), Strube (1971), Prosperetti (1977), Plesset and Prosperetti 
(1977), Brenner, Lohse, and Dupont (1995), Hilgenfeldt, Lohse, and Bren- 
ner (1996)), (ii) diffusional instabilities (Crum and Cordry (1994), Lofstedt, 
Barber, and Putterman (1993), Brenner et al. (1996a), Hilgcnfeldt, Lohse, 
and Brenner (1996)), and (iii) chemical instabilities (Lohse et al. (1996)). 
Hilgenfeldt, Lohse, and Brenner (1996) consider only the instabilities (i) and 
(ii). They could quantitatively understand the experimental results (Barber 
et al. (1994), Barber et al. (1995), Lofstedt et al. (1995)) for pure argon bub- 
bles, but not for gas mixtures with molecular constituents as e.g. air. For 
air the theoretical phase diagrams look very similar to those of argon, but 
experimentally (Barber et al. (1994), Barber et al. (1995), Lofstedt, Barber, 
and Putterman (1993)) stable SL is found for about 100 times larger gas 
concentration, as already mentioned above. In this paper we elaborate our 
idea (Lohse et al. (1996)) that considering also chemical instabilities resolves 
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this mystery. 



2 Phase Diagrams for Sonoluminescence 

To obtain phase diagrams of the bubble dynamics, one would have to solve the 
full three dimensional gas dynamical PDEs inside the bubble, coupled to the 
gas and fluid dynamics outside. This is numerically intractable; first, because 
the parameter space to be examined is huge, second, because a solution has 
to be computed for many millions of forcing cycles to cover the slow time 
scales of diffusional processes. 

Consequently, we do need approximations to continue. The most success- 
ful approximation of bubble dynamics is the Rayleigh-Plesset (RP) equation 
for the radius R(t) of the bubble as a function of time. This ODE was first 
written down by Lord Rayleigh (1917), while employed by the Royal Navy 
to investigate cavitation. Modern improvements have been made by Plcs- 
sct (1949), Lauterborn (1976), Plesset and Prosperetti (1977), Prosperetti 
(1977), and others. Recently, Fyrillas and Szeri (1994) and Lofstedt et al. 
(1995) were able to incorporate mass diffusion in the RP approach. This 
progress allowed us to study diffusional stability (Brenner et al. (1996a), 
Hilgcnfcldt, Lohse, and Brenner (1996)). 

The RP equation reads 

RR + |i? 2 - — (p(R, t) - P(t) - P ) 

2 Pw 

R d Ns R 2cr . . 

+ ^Jt^- P ^-^R-^R W 

with the forcing pressure field 

P(t) = Pa cos wt . (2) 

Parameter values for an argon bubble in water are the surface tension a = 
0.073kg/s 2 , water viscosity v = 10 _6 m 2 /s, density p w = 1000kg/m 3 , and 
speed of sound c w = 1481m/s. The driving frequency of the acoustic field is 
u/2ir = 26.5kHz and the external pressure P = 1 atm. We assume that the 
pressure inside the bubble varies according to the van der Waals law 

«"i=(4)^)*' (3) 

where h = i?o/8.86 is the hard core van der Waals radius. The effective 
polytropic exponent is taken to be k = 1 (Plesset and Prosperetti (1977)) for 
the small bubbles employed here. 

In figure 1 we plot R(t) resulting from (1), together with the forcing 
pressure P(t). The radius shows a characteristic collapse which can be very 
violent for large forcing P a . It is right at this collapse when the light pulse is 



4 Michael P. Brenner, SaschaHilgenfeldt, and DetlefLohse 




1 2 

time/T 



Fig. 1. (a) Forcing pressure P(t) = P a cosujt, P a — 1.3 atm for two cycles and 
the corresponding (b) bubble radius R(t) and (c) bubble ambient radius Ro(t). 
The bubble is near an equilibrium state. The gas concentration is c^, r / c o = 0.004. 
The ambient radius shown here results from the numerical solution of the convec- 
tion-diffusion PDEs, for details we refer to Hilgenfeldt, Lohse, and Brenner (1996). 



emitted. Barber and Puttcrman (1991) and Lofstcdt, Barber, and Puttcrman 
(1993) have shown that the experimental curves of the radius look very similar 
to that of figure 1 and that (1) can be used to fit the dynamics of the radius. 
There may be small quantitative discrepancies with the real dynamics, but 
they are clearly in the range of the experimental error. 

First, we address the question of shape instabilities , i.e., non-spherical de- 
formations of the gas-water interface. There are basically two different types 
of shape instabilities, distinguished by the time scale over which they act. 
Parametric instabilities act over many bubble oscillation periods, and are 
the result of accumulation of small perturbations. Rayleigh- Taylor instabil- 
ities occur when a light fluid is accelerated into a heavy fluid. Here, they 
happen during a short time directly after the bubble collapse. Although the 
process is very rapid (~ 10~ 9 s), detailed calculations (Brenner, Lohse, and 
Dupont (1995), Hilgenfeldt, Lohse, and Brenner (1996)) show that at high 
forcing, a perturbation of a few A can grow to be as large as the bubble 
during this time. The differing time scales over which these instabilities act 
have consequences for the final state of the bubble: after a Rayleigh- Taylor 
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instability, bubbles are typically destroyed, whereas after a parametric insta- 
bility, the bubble can remain trapped in the acoustic field. Figure 2 shows 
the borderline between stability and instability in the ambient radius Ro vs 
forcing amplitude P a parameter space. The diagram is calculated for an ar- 
gon bubble in water (Hilgenfeldt, Lohse, and Brenner (1996)). Also depicted 
in figure 2 is the threshold curve beyond which the bubble wall speed exceeds 
the speed of sound. We take this as a condition for the development of shock 
waves in the bubble. Following the common assumption that shocks are nec- 
essary for sonoluminescence to occur (Greenspan and Nadim (1993), Wu and 
Roberts (1993), Moss et al. (1994)), we thus get a further restriction for the 
parameter range of sonoluminescing bubbles. 




P/atm 

a 

Fig. 2. Parametric instability (solid), Rayleigh- Taylor instability (dashed), and 
shock condition (Mach=l, long dashed). SL is only possible in a small regime 
of large P a ~ 1.3 atm and small 7?o ~ 3/xm where the bubble is spherically sta- 
ble and the shock condition is fulfilled. Also shown are the diffusive equilibria for 
c^/co = 0.5, 0.02, and 0.002. Whenever the slope is positive, the diffusive equilib- 
rium is stable. At P a = 1.3atm diffusively stable SL bubbles are only possible for 
very low concentration c^ r /co ~ 0.002 — 0.004. For larger P a even lower concentra- 
tions are required. 

Let us now turn to diffusive stability. Measurements of the time be- 
tween successive light flashes show that the total mass of the bubble can 
remain constant to high accuracy (Gaitan et al. (1990), Barber and Putter- 



6 Michael P. Brenner, SaschaHilgenfeldt, and DetlefLohse 



man (1991), Barber et al. (1995), Lofstcdt et al. (1995)) for days, if the gas 
concentration is in an appropriate regime. This result contradicts classical 
notions about the dynamics of periodically forced bubbles: An unforced bub- 
ble of ambient radius i?n dissolves over a diffusive time scale, r ~ n , Pofi ° , 
(Epstein and Plcsset (1950)), where po is the ambient gas mass density in the 
bubble and D is the diffusion constant of the gas in the liquid. On the other 
hand, a strongly forced bubble grows by rectified diffusion, as first discovered 
by Blake (1949) and later explored by Eller (1969), Crum (1980), Crum and 
Cordry (1994). This is because at large bubble radius the gas pressure in the 
bubble is small, resulting in a mass flux into the bubble; conversely, when 
the bubble radius is small there is a strong mass outflux. Both processes can 
be observed in figure lc which shows the ambient radius Ro(t) corresponding 
to the actual bubble radius R(t) of figure lb. Since the diffusive time scale is 
much larger than the short time the bubble spends at small radii, gas cannot 
escape from the bubble during compression; the net effect is bubble growth. 
At a special value Rq of the ambient radius rectified diffusion and normal dif- 
fusion balance. The main result of Brenner et al. (1996a) was that contrary to 
above intuition the equilibrium points Rq can be stable. Hilgcnfcldt, Lohse, 
and Brenner (1996) showed that for argon bubbles the window of stability in 
c oo 7 c o is m quantitative agreement with experiment. 

Figure 3 shows the total phase diagram as a function of c^/co and P a / Pq. 
Three different states are possible: The label stable SL denotes the region of 
parameter space where the bubble can be diffusively stable and also sta- 
ble to shape oscillations. In the region no SL the bubble dissolves, while in 
the unstable SL domain the bubble is not in diffusive equilibrium but can 
nonetheless survive. That this is possible is quite nontrivial, and was first 
shown by experiments of Gaitan et al. (1990) and Barber et al. (1994). 

In the unstable SL phase the bubble is growing by rectified diffusion until 
it hits the parametric instability. The bubble becomes shape unstable and mi- 
crobubbles pinch off. If the remaining bubble is still large enough, the process 
begins again. Figure 4a shows the ambient radius Po an d the corresponding 
phase <j) s of the light emission relative to the forcing field as a function of 
time as calculated from our theoretical Rayleigh-Plessct approach (Hilgcn- 
feldt, Lohse, and Brenner (1996)). We chose the same three gas concentrations 
as in the experiment of Barber et al. (1995) which we present in Fig. 4b for 
comparison. Both in our calculations and in experiment we are in the un- 
stable SL regime for the two larger argon partial pressures (200 mmHg and 
50mmHg), whereas the lowest value of 3 mmHg is in the stable SL regime, 
again both in theory and experiment. Argon bubbles thus seem to be well 
understood. 

Another experimental result (now for a gas mixture) in the unstable SL 
regime is shown in Fig. 5. We observe the maximal radius P m ax and the light 
intensity growing on the diffusional time scale ~ Is until the bubble hits the 
parametric instability and the process starts over. Note that our interpreta- 
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Fig. 3. Phase diagram in the c^ r / c o vs P a /Po parameter space. For the SL phase 
we assumed that the light production mechanism works beyond the Mach=l curve. 
Bubbles at very large P a may be Rayleigh- Taylor unstable. 



tion of the unstable SL regime also accounts for the observed "dancing" of 
the bubble: the pinch-off of the microbubbles leads to recoils of the bubble. 
The dancing frequency should be the higher the larger c^ r is. 

Although the above results for pure argon bubbles are in excellent agree- 
ment with experiments, we still have above mentioned severe discrepancy for 
air bubbles and other gas mixtures that stable SL is observed for much larger 
concentrations than expected from our phase diagram Fig. 3. E.g., in air we 
have stable SL for c^ r / c o « 0.1 — 0.2, depending on the forcing pressure. This 
discrepancy was first pointed out in the qualitative arguments of Lofstedt et 
al. (1995) who found that air bubbles in the range where stability is observed 
should be unstable from a theoretical point of view, based on RP dynamics. 
The reason for this can be seen in the phase diagram above: if the diagram is 
repeated for air instead of argon, the stable equilibria still occur at very low 
concentrations, namely /c ~ 0.004. Based on this discrepancy, Lofstedt 
et al. hypothesize an "as yet unidentified mass ejection mechanism." 



8 Michael P. Brenner, SaschaHilgenfeldt, and DetlefLohse 



a) 





U 


CO 

3. 


-2 


-e- 


-4 




5.0 


B 




zL 


3.0 








1.0 



b) 




0.0 0.2 0.4 0.6 0.8 

t/s 

»■=-( -> 1 ! 1 , I 



- 



\ i 
U 

1 



i / 



L iff »»#* 




■1.44 



0.-1 



0.6 



Tim; ( s ) 



3mm Af & 

50mm Ar > 
— - 200mm Ar 

o 150mm Air 



—1— 

OS 



Fig. 4. (a) The phase of light emission 4> s (t) (upper) and the corresponding am- 
bient radius Ro{t) (lower) for P a = 1.3 atm for three different gas concentrations 
c^/co = 0.004, c^ r /co = 0.066, and c* r /co = 0.26, corresponding to a gas pres- 
sure of 3 mmHg, 50 mmHg, and 200 mmHg, respectively. The strength of the mi- 
cro-bubble pinch-off at Ro = 5/im, i.e. the decrease of the ambient radius, is picked 
randomly. The concentration values are chosen as in the experiment of Barber et 
al. (1995) which is shown in (b), copied from Fig. 4 of that reference. That fig- 
ure also shows the relative phase of light emission for air bubbles: Stable SL is 
achieved at much higher concentration c^> r / c o = 0.2, corresponding to a gas pres- 
sure of 150 mmHg. It is this discrepancy between argon (or other inert gases) and 
air which we resolve in this paper. 



Why air bubbles in water glow so easily 9 




Fig. 5. This figure is reproduced from figure 6 of Lofstedt et al. (1995). It shows 
the SL intensity and bubble radius as a function of time for a nitrogen bubble with 
5% argon. The drive level Pa. is near the upper threshold of SL. The gas is dissolved 
at 150mmHg, i.e. c™ lx /c ~ 20%. The bubble radius exhibits the characteristic 
pattern of growth and microbubble pinch-off of unstable SL. 

3 A Chemical Resolution of the Paradox 

In joint work with Todd F. Dupont and Blaine Johnston of the University 
of Chicago (Lohse et al. (1996)), we have pointed out that this paradox can 
be resolved by considering chemical instabilities. Namely, it is well known 
that the temperature in the bubble becomes very high after the collapse. 
Fits of the spectra of emitted radiation to a black body law give effective 
temperatures of up to 25000K (Hiller, Putterman, and Barber (1992)). This 
temperature exceeds the dissociation temperature for nitrogen gas (~ 9000K) 
as well as for all other molecular constituents of air. The nitrogen, oxygen, 
and hydrogen radicals will recombine to finally form NO, NO2, and/or NH3. 
All of these gases are highly soluble in water, forming nitric and nitrous acid 
and NH4OH. The consequence is that the bubble is finally depleted from the 
molecular constituents of air. 

Even beyond "burning" of all initial nitrogen in the bubble we expect 
that whenever the bubble is large it will still suck gas from the water into the 
bubble which will also be "burnt" during the next compression cycle. The 
bubble thus constitutes a reaction chamber for the dissolved molecular gas. 
The reaction products should be detectable. Lohse et al. (1996) estimated a 
production rate of ~ 10 11 mol/s per container volume (~ 11). Let us assume 
that the main product is NH4OH. Starting with the optimal pH=7 we should 
thus get a doubling of OH - ions within about three hours and a pH change up 
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to pH=8 (i.e., c OH - = 10~ 6 mol/l) in about a day of running the experiment. 
Many other ways of detection rather than simply measuring the pH seem 
possible. An ion chromatographcr will give a detailed answer on what ions 
are produced. 

For the present paper we only need the assumption that even when air 
was initially dissolved in water, a sonoluminescing bubble essentially consists 
of inert gases as these are the only gases which do not react with water at the 
high temperatures achieved in single bubble SL; they are simply recollected 
by the expanding bubble as seen in Fig. lc. 

The dissociation hypothesis suggests that for gas mixtures the relevant 
concentration quantity is still the argon concentration c^/cq which now is 



where q is the percentage of argon in the mixture. For air we have q = 0.01 = 
1%. Here we have neglected the differences in solubility Co of the different 
gases. This difference enters in two counteracting ways: first, (4) must be 
corrected by the ratio Cq t /Cq 2 of the solubilities between argon and nitrogen, 
which is about three. But second, because of the different solubilities the 
percentage of argon in dissolved air will not be the usual 1%, but will be 
larger. Possibly, this percentage even depends on the degree of degassing of 
the water. These two effects will probably roughly compensate each other 
and we decided to disregard them for the time being. 

In the following section we discuss how the dissociation hypothesis can 
explain the experimental UCLA results on gas mixtures (Hiller et al. (1994), 
Barber et al. (1994), Barber et al. (1995), Lofstedt et al. (1995)). We point 
out that our detailed unpublished investigations of mixtures of gases demon- 
strate that it is impossible to understand these discrepancies within diffusive 
theories alone, i.e., by disregarding chemical reactions. 

4 Comparison with the UCLA Experiments 

Where to expect stable SL in air bubbles? From our phase diagram Fig. 3 
we know that for P a = 1.3 atm stable SL is found in a small window of 
c oo/ c o ~ 0-3 — 0.4% for pure argon gas. Equation (4) thus gives c^ r / c o ~ 
30 — 40% for air dissolved in water. The experimental values seem to be 
slightly smaller, but clearly within the error range of both our assumptions 
and the measurements. 

The transition towards SL with increasing forcing pressure P a is shown in 
Fig. 6 for both argon and air bubbles. For pure argon bubbles the transition 
to SL is very smooth. For air, however, one can observe a breakdown of the 
bubble radius at about 1.1 atm, signaling that the dissociation threshold of 
N2 is achieved. Before the transition the bubble is filled with a mixture of 
nitrogen, oxygen, and argon, and the ambient radius is determined by the 
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Fig. 6. Transition towards the SL regime for argon bubbles (a) and for air bubbles 
(b). Only for air bubbles a breakdown in the radius is seen near the onset of SL, 
signaling the threshold for nitrogen dissociation. Figure (a) is a reproduction from 
figure 4 from Hiller et al. (1994), the total gas concentration is about 150mmHg. 
This means c^ r /c = 0.2 in the unstable SL regime in agreement with the obser- 
vations. Figure (b) is a reproduction from figure 2 of Barber et al. (1994). The gas 
saturation is about 10%, corresponding to c^ r /co = 0.01 • 0.1 = 0.001. According 
to Fig. 3 we have stable SL around P a ~ 1.4 atm which again is in rough agreement 
with experiment. 
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combination of all three gases. After the dissociation threshold, only argon is 
left in the bubble. 

We now also see why stable SL is so easily obtained in air bubbles: The 
window of stability is larger compared to argon bubbles and the water has 
to be degassed much less. But even more friendly conditions are possible: 
Our theory suggests that it should be possible to obtain stable SL without 
degassing, namely by choosing the percentage q of argon so that the window 
of stable SL is around c™ IX /co ~ 100%. With q ~ 0.35% we obtain as window 
of stable SL c™ ix /c - 85 - 115%. 




Fig. 7. This figure is reproduced from figure 11 of Lofstedt et al. (1995). It shows 
the transition to SL for a bubble filled with an initial 0.1% xenon in nitrogen gas 
mixture at a partial pressure of 150 mmHg. According to our dissociation hypothesis 
this corresponds to c^°/co = 0.0002 for a pure sonoluminescing xenon bubble. From 
our phase diagram Fig. 3 which is (with tiny corrections) also valid for xenon we 
conclude that bubbles driven at P a = 1.3 atm dissolve for these low concentrations, 
whereas bubbles at P a = 1.4 atm show stable SL, just as seen in this experimental 
figure. 

For even lower concentration q < 0.3% (at P a = 1.3 atm) no stable SL 
regime is left for c™ IX /co < 100%. The experiment of Lofstedt ct al. (1995) 
with q — 0.1% xenon in nitrogen at c™ IX /co = 0.2 shown in Fig. 7 relates to 
exactly that situation. Lofstedt et al. (1995) find a range of forcing pressures 
P a « 1.3 — 1.4 atm where stable bubbles cannot be seeded. Above 1.4 atm, 
stable SL exists, and below 1.3 atm there are non-sonoluminescing bubbles. 
The reason for this "gap" is that at these high forcing pressures all nitrogen is 
dissociated and the bubble only contains the inert gas xenon. Then, according 
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to our phase diagram Fig. 3 (which also holds for xenon), at P a = 1.3 atm 
the xenon concentration c* c = q ■ c™ IX /c = 0.001 • 0.2 = 0.0002 is in the no 
SL regime and bubbles dissolve. But at higher forcing pressure P a ~ 1.4 atm 
we can again enter the stable SL regime for that concentration, see Fig. 3: a 
pure xenon bubble can emit light and acts as a chemical reaction chamber, 
transferring the dissolved nitrogen to its high temperature reaction products. 




Noble gas in N 2 (%) 



Fig. 8. This figure is reproduced from figure 1 of Hiller et al. (1994). It shows 
the SL intensity from a SL bubble in water as a function of the percentage (mole 
fraction) of noble gas mixed with nitrogen. The gas mixture was dissolved in water 
at a pressure head of 150mmHg, i.e., c™ lx /co ~ 0.2. 

The transition from the no SL regime in Fig. 3 to the SL regime can 
also be seen in Fig. 8, where the SL light intensity is plotted as a func- 
tion of q for fixed P a (we assume P a = 1.3 atm) and fixed c™ lx /co ~ 0.20. 
According to our theory the transition for c^ r / c o = 0.003 should occur at 
q = (c^ r / c o)/( c ™ lx / c o) = 1.5% in pretty good agreement with Fig. 8 where 
we indeed observe that strong SL is "switched on" at about that concentra- 
tion. Our theory also predicts that near the switch on we always have stable 
SL, whereas for larger q unstable SL develops. Here we expect unstable SL 
for q > 2%, a prediction which should be verified. 

An example for unstable SL was already shown in Fig. 5. Indeed, for that 
figure we have q = 5%. Thus, from (4) we have c^/c = 0.05 • 0.2 = 0.01 
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and according to our phase diagram Fig. 3 and equation (4) we are well in 
the unstable SL regime, just in agreement with the observations. 

Next, we consider the case of a pure nitrogen bubble. Below the dissocia- 
tion threshold it is "jiggling" , i.e. growing by rectified diffusion and shedding 
microbubbles. Above the threshold, the nitrogen is taken away, so the bubble 
will necessarily dissolve. Indeed, Hiller ct al. (1994) show that pure nitrogen 
bubbles dissolve in a matter of seconds. In figure 5 of their work it is demon- 
strated that pure nitrogen bubbles show "periodic" oscillations in their light 
intensity. We do not understand what sets the oscillation period, though we 
speculate that it is connected to the dissociation of nitrogen in the bubble, 
as well as the unstable "jiggling" of the bubble. Perhaps by understanding 
the origin of these oscillations one could deduce the temperature inside the 
bubble. 

Another experiment supports our dissociation hypothesis. Hiller and Put- 
terman (1995) analyze SL in H2 and D2 gas bubbles, both in normal and in 
heavy water. One would expect that the SL intensity curves would group ac- 
cording to the gas content as it is the gas dynamics inside the bubble which 
determines the strength of the light emission. However the four (H2 in H2O, 
H2 in D 2 0, D 2 in H 2 0, and D 2 in D 2 0) experiments group according to 
the surrounding liquids. This suggests the following scenario: Both the gas 
and the liquid vapor in the bubble dissociate during the (hot) compression 
phase and recombinc later on during expansion. As much more liquid than 
gas is around, recombination will yield gases primarily composed of the con- 
stituents of the liquid vapor. E.g., in the experiment using H2 in D2O, the 
H2 will shortly be replaced by D 2 formed from the heavy water, whereas the 
amount of H2O formed from H2 will be negligible. Thus the density of the 
bubbles in D2O will eventually be larger by a factor of two compared to the 
density of bubbles in H 2 0. Following our recently proposed laser theory of 
SL (with R. Rosalcs, (Brenner et al. (1996b))), it is this density difference 
which may be the origin of the considerably different SL intensities observed 
in normal and heavy water by Hiller and Putterman (1995). 

There is one final experimental figure which we want to discuss here. It is 
Fig. 9, showing the maximal radius i? ma x(i) = max{R(t')\t < t' <t + T} and 
the intensity as a function of time when the forcing pressure P a is suddenly 
strongly increased. On this increase, the bubble is pushed in the unstable SL 
domain and starts to grow rapidly thanks to rectified diffusion. But why is 
the growth rate of R ma , x so wiggly whereas the growth rate of the ambient 
radius does not show these wiggles, as seen from Fig. 4? 

The answer is that i? max itself shows wiggles as a function of the ambi- 
ent radius, as observed in Brenner et al. (1996a) and analyzed in detail by 
Grossmann et al. (1996). In figure 10 we present R ma , x (Ro), whose wiggles 
are probed in figure 9 by the growing R n . The averaged scaling law is given 
by -Rmax oc Rq^ 5 (Grossmann et al. (1996)). 

In diffusionally stable regimes these wiggles may lead to multiple stable 
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Fig. 9. This figure is reproduced from figure 4 of Barber et al. (1994). It shows the 
dynamic response of the radius and the intensity of SL to a sudden change in the 
drive level P a . The bubble is boosted into the unstable SL regime where it finally 
becomes shape unstable and bursts. 

equilibria, as suggested in Brenner et al. (1996a) and possibly observed by 
Crum and Cordry (1994). Grossmann et al. (1996) pointed out that the origin 
of the wiggles and the resulting multiple equilibria are resonances in the RP 
dynamics. 

5 Conclusions 

This paper has summarized our current understanding of the bubble dynam- 
ics leading to single bubble sonoluminescence. Our central theoretical result 
is the phase diagram of SL figure 3 from Hilgenfeldt, Lohse, and Brenner 
(1996). It can only be calculated in the Rayleigh-Plcsset SL bubble approach 
including both shape oscillations and diffusion. The phase diagram allows to 
quantify the differences between argon and air bubbles found by Lofstedt et 
al. (1995), guiding them to suggest a "yet unidentified mass ejection mech- 
anism" for air bubbles. We now believe that this mechanism is chemical, 
namely the dissociation of gas molecules in highly forced gas bubbles and 
their subsequent reactions. As pointed out in this paper, with this "nitrogen- 
dissociation hypothesis" many experimental results can be explained. To fur- 
ther verify this hypothesis, experiments should search for the various reaction 
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Fig. 10. The maximal bubble radius i? m ax as a function of the ambient radius Rq 
for two forcing pressures P a = 1.1 atm and 1.2atm. It shows characteristic wiggles. 
The dashed lines depict the -R max oc R^ 5 scaling law. 



products. 

One issue that we have barely touched in this paper is the mechanism 
for the light emission. The traditional theory is that shocks are emitted from 
the bubble wall, and heat up the gas in the bubble to such an extent that 
light emission is induced. Experimentally there are strong dependences of the 
light intensity on liquid, gas and forcing pressure. These cannot be explained 
within the context of simple bubble dynamics presented here, and also do 
not seem to be understandable within the shock theory. The shock theory 
relies on a singular solution to the hydrodynamic equations; singularities are 
typically of universal character and very insensitive to external conditions. 
However, the fact that experiments on single bubble SL demonstrate such sen- 
sitivity to parameters contradicts this principle. The most dramatic example 
of this is the initial observation of single bubble SL by Gaitan et al. (1990): it 
caused great excitement because the light intensity is so much higher than in 
multibubble SL. However, the shock theory should apply equally well to both 
cases. These ideas suggest that the shock theory has intrinsic shortcomings. 
To overcome these problems Brenner et al. (1996b) recently proposed a "hy- 
drodynamic laser" theory for SL bubbles. It stipulates that the bubble can 
accumulate acoustic energy inside the bubble which acts as a kind of acoustic 
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resonator. 

We hope that our present and ongoing work will help to answer one of 
the essential remaining questions on SL: How hot does the gas become in 
the center of the bubble? Our hope is that the presence of intermediate or 
final chemical reaction products can be thought of as a thermometer. - The 
perspectives are appealing: A tiny air bubble in water is used as a chemical 
reaction chamber and a micro-laboratory for high temperature chemistry! 
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